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Abstract 

In this paper we describe two new computational operators, called complex en- 
tropic form (CEF) and generalized complex entropic form (GEF), for pattern char- 
acterization of spatially extended systems. Besides of being a measure of regularity, 
both operators permit to quantify the degree of phase disorder associated with a 
given gradient field. An application of CEF and GEF to the analysis of the gradient 
pattern dynamics of a logistic Coupled Map Lattice is presented. Simulations using 
a Gaussian and random initial condition, provide interesting insights on the the 
system gradual transition from order /symmetry to disorder /randomness. 
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1 Introduction 

Nonlinear spatially extended dynamical systems yield complex amplitude pat- 
terns which arise from the coupled dynamics of their different regions. Exper- 
iments in a variety of areas have shown spatio-temporal complexity, notably 
in fluid flows, diffusive-reactive chemical systems, optical electronics and laser 
physics. As reported by many authors (e.g. [1]), these systems cover a wide 
range of scales, from millimeters in optical electronics experiments to thou- 
sands of kilometers in natural systems. In all this cases it is important to 
quantify the degree of local complexity in order to characterize the spatial 
patterns and also study their time evolution. Consequently, one of the main 

1 Corresponding author: F.M.R., fernando@lac.inpe.br 
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tasks in computational physics today is to define suitable measures to char- 
acterize the dynamics of such complex systems. As stressed by Chate [2], the 
current understanding of spatio-temporal disorder in extended systems is still 
very limited due to the lack of quantitative, meaningful, and sufficiently uni- 
versal patterns and regimes characterizations. 

Recently, the analysis of the gradient field of the amplitude envelope of dynam- 
ical extended systems has been shown to be a useful tool for understanding 
complex regimes as intermittency and localized turbulence [3]. However, the 
classical measures of complexity usually discharge the directional information 
contained in a vector field. In this paper we describe two new computational 
operators, called complex entropic form (CEF) and generalized complex en- 
tropic form (GEF), for pattern characterization of spatially extended systems. 
Besides of being a measure of regularity, both operators permit to quantify 
the degree of phase disorder associated with a given gradient field. Phase dis- 
order plays an important role in the analysis of pattern changing frequency in 
the amplitude domain [4]. For illustration purposes, we apply CEF and GEF 
to the analysis of the gradient pattern dynamics of a logistic Coupled Map 
Lattice. 



2 Complex Entropic Form (CEF) 

If we consider spatially extended systems in two dimensions (x,y), their energy 
amplitude distribution is represented by the envelope A(x,y), which can be 
approximated by a matrix of amplitudes A = {a^i}, with K x L pixels. Note 
that a dynamical sequence of matrices can be related to a temporal evolution 
of an envelope A(x,y,t). From the definition of A, it is possible to represent the 
gradient field of the amplitude envelope by G = V(A) = {zk,i}, where Zk,i is a 
complex number with Re(z k> i) = a k ,i+i - a k ,i-i and Im(z k j) = a k+ i,i - a k -i,i- 

Given the gradient field matrix V(A), we now define the Complex Entropic 
Form (CEF) operator as: 

S C (G) = -EW = - E W H %r) - i E tjt^m (i) 

where \z\ = J2k,i\ z k,i\, \z k ,i\ and 0*^ are respectively the modulus and the 
argument of z k ,i and i = T- 

From the above definition, we immediately verify that the Re(S c ) corresponds 
to the classical Shannon's entropy measure S of the moduli \zk,i\ and that 
Im(S c ) represents the weighted average phase of the gradient field. Conse- 
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quently, a random generated pattern whose gradient field displays no domi- 
nant direction, will have Im(S c ) ~ and Re(Sc) close to its maximum value, 
Re(S c ) = log(M), where M = K x L is the number of elements in the gradient 
matrix. On the other hand, intermittencies and symmetry breaking patterns 
lead to, respectively, a decrease on the value of Re(S c ) and to a non null 
Im(S c ). 



3 Generalized Complex Entropic Form (GEF) 

A more general entropic form that also incorporates the phase information 
of the gradient field can be obtained from a generalization of the concept of 
degeneracy W, given by the multinomial coefficient formula, and normally 
used to deduce the expression of Shannon's entropy of positive scalar fields, 
such as images [5]. Considering the gradient matrix G defined in the previous 
section, W(z lt i, . . . , zk,l) ma Y be generalized as follows: 

w{ *» «" )m nJ.% K *) ' (2) 

where z = J2k,i z k,i- Using Stirling's approximation, we immediately have 

z -HnW ^ S z = -J2 — H-) • (3) 

k,i z z 

This expression, that we hereafter call Generalized Complex Entropic Form 
(GEF), displays interesting properties [6], several of which it shares with 
Shannon's entropy. In particular, we can easily verify that S z is invariant 
under rotation and scaling of the vector field, and that Im(S z ) = and 
Rz{S z ) = S(\zk t i\), when all z^j, have the same direction (in other words, 
no information is conveyed by the phases when they are all the same). 



4 Application 

In order to illustrate the performance of the computational operators described 
in the previous sections, we applied CEF and GEF to the characterization of 
amplitude pattern formation (fragmentation, symmetry breaking and phase 
disorder) of a logistic extended Coupled Map Lattice (CML). Additional re- 
sults of the performance of these operators in different contexts (osmosedi- 
mentation process [7], wavelet spectra of temporal series [8], Swift-Hohenberg 
dynamics [9], porous silicon surface patterns [10]) are reported elsewhere. 
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CML is one of the most important systems that have been introduced to study 
the dynamics of spatio-temporal complexity in the amplitude domain [11]. The 
explicit form of a global CML is given by 

M 

x n+1 = (1 - e)f(x n (k)) + -E f(x n (l)) , (4) 

where n refers to a discrete time step and M = K = L is the lattice size. 
Usually, the function f(x) is chosen to be the well known dissipative chaotic 
logistic map, f(x) = 1 — px 2 . 

From Eq. (4), considering the matrix of amplitudes A = {a^}, we have 

<f = (1 " eK, + + + f(aU,i) + , (5) 

where e is the coupling constant and f(a) = pa(l — a), with p being a real 
constant (the chaotic control parameter). 

We solved Eq. (5), for e = 0.5, p = 4.0 (chaotic regime), and an aspect ratio 
(ratio between largest and smallest system characteristic scales) of T = 400. 
We started the simulations from two different initial conditions: (i) a Gaus- 
sian, symmetric initial distribution, and (ii) a random distribution. Figure 1 
(a,b,c,d) shows four frames (n = 5, 20, 35, 50) of the evolution of the gradient 
field and the corresponding intensity contours plots for the Gaussian initial 
condition. Figures 2 and 3 depicts respectively the variation of Re(S z ) and 
Re(S c ), for test cases (i) and (ii). 

While the results of the random test case, as expected, do not display any 
particular characteristic, from the evolution of Re(S z ) for the Gaussian initial 
distribution (Figure 2) we identify three well defined periods, corresponding 
to different phase disorder pattern regimes. For n < 5, the value of Re(S z ) is 
very high (but fast decreasing), since the gradient field still presents a high 
degree of symmetry. For 5 < n < 20, as the Gaussian perturbation reaches the 
boundaries of the system, the symmetries of the initial condition are almost 
completely lost. Moreover, during this period, the system have not yet acquired 
the symmetries, in a statistical sense, that characterizes a random field and, 
thus, Re(S z )c < Re(S z )n. Finally, for n > 20, the gradient field becomes 
increasingly random. However, since the system still carries a memory of the 
original symmetry, the value of Re(S z )c is higher than Re(S z )n (see [12] for 
a similar analysis, concerning the loss and recovery of symmetries of Navier- 
Stokes equation for increasing values of Reynolds number). The evolution of 
Re(S c ) (Figure 3), although with less detail, also shows the system gradual 
transition from order/symmetry (lower Re(S c ) value) to disorder /randomness 
(higher Re(S c ) value). 
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5 Concluding Remarks 



In this paper we introduced two new computational operators, called com- 
plex entropic form (CEF) and generalized complex entropic form (GEF), for 
pattern characterization of spatially extended systems. The real part of CEF 
corresponds to the classical Shannon's entropy measure, while its imaginary 
part represents the weighted average phase of the gradient field. The definition 
of GEF is obtained from a generalization of the concept of degeneracy. Besides 
of being a measure of regularity or smoothness of the embedding amplitude 
envelope, both operators provide information on the degree of phase disorder 
of the corresponding gradient field. 

For illustration purposes, we applied CEF and GEF to the analysis of the gra- 
dient pattern dynamics of a logistic Coupled Map Lattice. Simulations using 
a Gaussian and random initial condition, provided interesting insights on the 
the system gradual transition from order/symmetry to disorder /randomness. 
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Figure Captions 



Fig. 1 Four frames (n = 5, 20, 35, 50) showing the evolution of the gradient 
field and intensity contours plots for the solution of logistic Coupled Map with 
Gaussian initial condition, T = 400, e = 0.5 and p = 4.0. 

Fig. 2 Evolution of Re(S z ) as a function of the time index n, for Gaussian 
and random initial conditions. 

Fig. 3 Evolution of Re(S c ) as a function of the time index n, for Gaussian 
and random initial conditions. 
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